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Abstract 

We compute the dilatation generator in the su(2) sector of planar Af = 4 
super Yang-Mills theory at four-loops. We use the known world-sheet 
scattering matrix to constrain the structure of the generator. The re- 
maining few coefficients can be computed directly from Feynman dia- 
grams. This allows us to confirm previous conjectures for the leading 
contribution to the dressing phase which is proportional to C(3)- 
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1 Introduction and Overview 

The means available for analyzing the AdS / CFT correspondence improved dramatically 
with the discovery of perturbative integrability of the gauge theory dilatation operator 
[IHS] and that of classical integrability of the world sheet sigma model [1]. Furthermore, 
there are arguments [5J on the string theory side of the correspondence that a infinite 
family of BRST invariant, non-local currents exists at all orders in the inverse 't Hooft 
coupling expansion suggesting that integrability persists in the quantum theory. In 
the absence of a definitive and constructive proof of all-order integrability, one may 
nonetheless assume it and study its consequences. 

The fundamental quantity in an integrable (discrete or continuous) theory defined 
on infinitely extended space-like slices is the scattering matrix of excitations. The S- 
matrix is constrained by the symmetries of the theory; integrability further requires that 
no particle production occurs in the scattering process and that the n —>■ n scattering 
process is realized by repeated 2 — > 2 scattering events. A necessary requirement is that 
the two-particle S-matrix obeys the Yang-Baxter equation. 

For the AdS/CFT correspondence the relevant two-particle scattering matrix was 
introduced in [6]; it turns out that the global symmetries - a centrally extended form 
of psu(2|2)^ determine it up to an overall phase [7]. The Yang-Baxter equation holds 
automatically in this case. Although initially the S-matrix was determined in the gauge 
theory framework it was later shown that the tensor structure agreed with perturbative 
calculations in the gauge-fixed world-sheet theory [8] and that it is consistent with the 
Zamolodchikov-Faddeev algebra for the string sigma- model p]. 

In relativistic quantum field theories the analogous "dressing factor" is determined 
by crossing symmetry, information on the spectrum of bound states and perhaps pertur- 
bative calculations. For the AdS/CFT correspondence both the world sheet and gauge 
theory integrable systems do not exhibit Lorentz invariance. While on the gauge theory 
side there is little reason to require an analog of crossing symmetry, on the string theory 
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side two-dimensional Lorentz invariance is only spontaneously broken. As such, one may 
expect that some form of crossing symmetry survives this breaking. 

A crossing-like equation was constructed in [TU] and shown in [TT] to hold for the 
known leading [T2j and next-to- leading terms [13]. An all-orders solution at a strong 
coupling expansion was proposed in [H]. 

An unfortunate feature of this solution is that it is an asymptotic series and thus, 
without additional information, cannot be directly used to define the dressing phase 
everywhere in the coupling constant space. In [IS] an analytic continuation scheme was 
described which allowed a guess for the weak-coupling expansion of the dressing phase 
whose contribution to anomalous dimensions starts at four-loop order where it predicts a 
transcendental contribution proportional to C(3)- This prediction remarkably agrees with 
the direct calculation of the four-loop cusp anomalous dimension [T6l[I7]. Subsequently 
the expansions at weak and strong coupling were shown to be fully consistent [IH] and 
an integral expression for the phase at finite coupling was proposed in [TH] 

In fact the above agreement is slightly surprising: The four-gluon scattering am- 
plitude of [16] is related to the infinite-spin limit of twist-two anomalous dimensions. 
Conversely, the analysis of [15j strictly applies to local operators of twist three or higher. 
Due to the asymptotic nature of the higher-loop Bethe equations the twist-two anoma- 
lous dimension can only be predicted reliably up to three loops, see [SU] for further recent 
developments. The agreement thus implies that the cusp anomalous dimension is univer- 
sal for operators of all twists. In other words, the limiting procedure described in [2T1[T5] 
does not suffer from potential order-of-limits ambiguities. 

As remarkable as it is, this agreement also presents a puzzle: The universality of 
the dressing phase implies that all anomalous dimensions of A/" = 4 SYM have, at four- 
loop order, a transcendental contribution proportional to C(3)- While this is not at all 
surprising for non-compact subsectors of A/" = 4 SYM in the large spin limit, it does 
seem surprising for finite spins and for compact sectors. Indeed, in the infinite spin limit 
the RG flow mixes an infinite number of operators allowing transcendental numbers to 
appear even if they are absent at the level of the anomalous dimension matrix. In the 
latter cases however, the RG flow mixes only finitely many operators and thus precludes 
the appearance of transcendental numbers. Consequently, for the conjectured dressing 
phase to be correct, ({3) must appear at the level of the anomalous dimension matrix 
elements. 

Loop integrals may be interpreted - in a first quantized language - as a sum over 
infinitely many intermediate states producing an analogy with the large-spin s[(2) sec- 
tor operator mixing. From this standpoint, one is entitled to expect the appearance 
of transcendental numbers at some sufficiently high loop order in any sector. One of 
the building blocks of the calculation of the renormalization factors of scalar composite 
operators is the one- loop scalar bubble diagram. It turns out that, in dimensional reg- 
ularization, its e-expansion contains ({3) at 0(e^); consequently, if this bubble is part 
of a larger diagram and the other momentum integrals yield a third-order pole in the e 
expansion, ^(3) may appear in the residue of a first order pole and thus may contribute 
to some entry of the anomalous dimension matrix. Counting the required number of 
inverse powers of the dimensional regulator we immediately reach the conclusion that 
this mechanism may function first at the four-loop order. 
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In this paper we shall compute the four-loop dilatation operator in the 3u(2) sector 
and show that the expectations outlined above are indeed realized. We shall begin in ^ 
with a review of the constraints imposed by su(2) symmetry and Feynman diagrammat- 
ics. The unknown coefficients are parametrized in terms of the first nontrivial coefficient 
of the dressing phase. However, unlike earlier discussions [22l|23] we shall not assume 
that this operator is part of an integrable Hamiltonian. Instead, we shall determine in 
§3] the unknown coefficients - and in particular the coefficient related to the dressing 
phase - by a direct calculation. The calculation is dramatically simplified by the obser- 
vation that the unknown coefficients may be associated to so-called maximal interactions 
(i.e. interactions that reshuffie the spins in a maximal way). §1] contains our conclusions. 
Some technical details as well as some momentum integrals useful for going beyond four 
loops are included in the appendices. 

2 Long-Range Heisenberg Hamiltonian 

A full-fledged fleld theory calculation of the complete four-loop planar dilatation genera- 
tor in A/" = 4 SYM is a difficult task whose completion clearly requires new, deep insight 
in higher-loop technology. The main complications are the extensive combinatorics and 
the intricate algebra of loop momenta inherent to gauge theories at higher perturbative 
orders. However, our primary goal is to compute the relevant coefficient for the dressing 
at this order. The dressing factor can be observed in all closed sectors of the model and 
we can conveniently restrict to the simplest one, the su(2) subsector, cf. [24J. It consists 
of local operators which are made from just two complex scalars, let us denote them 
by Z and 0, or, equivalently, spin up and spin down. Here the planar dilatation oper- 
ator turns into the Heisenberg XXX1/2 Hamiltonian [1] with perturbative long-ranged 
deformations [3] 



Determining this Hamiltonian at the fourth perturbative order would provide us with 
the leading piece of the dressing phase. 

The first few perturbative deformations of the Hamiltonian were obtained in [3]: 
This construction made use of the fact that the Hamiltonian is some linear combination 
of all interactions compatible with su(2) symmetry which can originate from Feynman 
diagrams. The coefficients of the interactions could in principle be computed from per- 
turbative field theory. However, such an elaborate calculation was avoided by matching 
the coefficients to make the spectrum of the Hamiltonian agree with some available data. 
Together with the further assumption of integrability, a proposal for the Hamiltonian at 
the third perturbative order could be made. The conjecture has since passed various 
tests [25l[2T] which prove that it is correct. 

Here we shall repeat the above procedure to constrain the fourth-order Hamiltonian 
as much as possible without making any unproven assumptions. The crucial new input 
that allows us to go to higher orders is the picture of asymptotic excitation states [6] and 
its scattering matrix [7]. In this picture, spin chain states are replaced by excitations 




(2.1) 
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above a ferromagnetic vacuum, the magnons. The ferromagnetic vacuum consists of a 
long chain of ahgned spins, say Z 

\Q) = \...ZZZ ...). (2.2) 

This state is protected by a half-BPS condition from receiving quantum corrections to 
its energy; the complete cancellation of corrections to two-point functions in field theory 
at two loops is demonstrated explicitly in [2S]- A single-magnon state has one of these 
spins flipped to 0, say at position k 

k 

I 

\k) = \...Z(j)Z ...). (2.3) 
Similarly, one can construct states with two or more magnons 

k I 

\k,l,..) = \. . . Z(t)Z . . . Z(t)Z . . . Z(f)Z ). (2.4) 

In the asymptotic coordinate space Bethe ansatz [6] the magnons are arranged into 
momentum eigenstates with an additional phase shift when two magnons move past 
each other. 

The excitation picture is highly constrained by its residual symmetry. It was shown 
in [7] that the form of the one and two-magnon states is almost completely determined. 
The only degrees of freedom are a finite redefinition of the coupling constant and the 
dressing phase. The possibility to redefine coupling constants by a finite amount is 
inherent to field theories. We can make a suitable choice and all other choices can be 
recovered from it by substitution. A general analysis [23] shows that the dressing phase 
starts to contribute at four loops with a single undetermined coefficient ^2,3- 

Note that these results are actually not based on the (unproven) assumption of higher- 
loop integrability: integrability or factorized scattering constrains the scattering of three 
or more particles. It also implies a constraint on the two-particle scattering matrix which, 
however, in this case is automatically satisfied [7]. 

We can now match the coefficients of the Hamiltonian to the zero-, one- and two- 
particle states. The analysis proceeds along the lines of [23] and the most general result 
is shown in Tab. [H The interaction symbols {a, b,c, . . .} represent a sequence of nearest- 
neighbor interactions Vp of spins at sites p and p + 1 summed homogeneously over the 
spin chain of length L 

L 

{a, 6, c, . . .} = ^ Vp+aVp+bVp+c ■■■ . (2.5) 
p=i 

As undetermined parameters it contains the coefficient /32,3 for the dressing phase as well 
as several irrelevant parameters e. The latter correspond to similarity transformations 
of the Hamiltonian which do not affect its spectrum. One may change their values by 
applying the similarity transformation Ti 1— ex'p{—iX)T-C exp[+iX) with the second and 
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^0 = +{} 

7^i = +2{}-2{l} 

^2 = -8{} + 12{1} - 2({1, 2} + {2, 1}) 

^3 = +60{} - 104{1} + 4{1, 3} + 24({1, 2} + {2, 1}) 

- Aie^il, 3, 2} + 4ze2{2, 1, 3} - 4({1, 2, 3} + {3, 2, 1}) 

7^4 



+ 



+ 



+ 
+ 
+ 
+ 



+ 
+ 
+ 



-560 - 4/52,3) {} 
+ 1072 + 12/?2,3 + 8e3a){l} 
-84 - 6/^2,3 -4e3a) {1,3} 
4{1,4} 

-302 - 4/32,3 -8e3a) ({1,2} + {2,1}) 

+4/32,3 + 4e3a + 2^63, - 4ze3d) {1, 3, 2} 

+4/32,3 + 4e3a - 2ze3c + 4ie3d) {2, 1, 3} 

4-2ze3c)({l,2,4} + {l,4,3}) 

4 + 2^630) ({1,3, 4} + {2, 1,4}) 

+96 + 4e3a)({l,2,3} + {3,2,l}) 

-12-2/32,3-4e3a){2,l,3,2} 

+ 18 + 4e3a)({l,3,2,4} + {2,l,4,3}) 

-8 - 2e3a - 22e3b) ({1, 2, 4, 3} + {1, 4, 3, 2}) 

-8 - 2e3a + 2ze3b) ({2, 1, 3, 4} + {3, 2, 1, 4}) 



-10({1,2,3,4} + {4,3,2,1}) 



Table 1: The four-loop Hamiltonian. The coefficient /32,3 is the leading coefficient of the dressing 
phase at weak coupling. We confirm the prediction /32,3 = 4(^(3) [15J as the principal result of 
this paper. The coefficients e correspond to similarity transformations and do not influence the 
spectrum. 
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third-order contributions to X given by 



A'2 = (5e2({l,2} + {2,1}), 

^3 = 25e3a({2, 1, 3} - {1, 3, 2}) + 5e3b ({1, 2, 3} + {3, 2, 1 
+ 5e3c{l,3} + 5e3d({l,2} + {2,l}). 



}) 



(2.6) 



It is worth pointing out that the Hamiltonian at fourth order can be fixed uniquely 
up to irrelevant terms. In other words, the scattering of three or more magnons is 
fixed by the scattering of two magnons. This feature is related to the su(2) symmetry 
of the interactions: Interactions at four loops act on at most five adjacent spins. Any 
elementary interaction among three or more magnons (and therefore at most two vacuum 
spins) is related to an interaction among at most two magnons (and at least three vacuum 
spins) by fiipping all five interacting spins. Starting at five loops this picture breaks down 
because interactions of six spins allow for elementary interactions of three magnons which 
leave no trace on the sector with two or fewer magnons. It turns out that our four-loop 
Hamiltonian in Tab. [U is integrable, i.e. it is of the form determined (but not displayed 
explicitly) in [23]. We have therefore proved four-loop integrability in the su(2) sector. 

The four-loop Hamiltonian in Tab. [T] is fixed to a large extent. To determine the 
dressing phase coefficient (32,3 it suffices to compute only a small number of its coefficients. 
We see that j32,3 couples, among others, to the very first and fifth but last interaction 
structure in Tab. [H The first structure does not redistribute the spins along the spin 
chain. There are exceedingly many planar Feynman diagrams which do not change 
flavor, for example those containing only interactions of gluons and scalars. Therefore, 
a direct computation of this coefficient seems particularly difficult. In contrast, the 
coefficients of the five last interactions can be computed relatively easily. They form a 
class of interactions which reshuffle the spins in a maximal way. At i loops, they contain 
i permutations of nearest neighbors, see Fig. [1] for a graphical representation of their 
induced permutations. This is the maximum reshuffling allowed by planar Feynman 
diagrams |3J and it will turn out to be generated by the quartic interactions of the scalars 
only. In other words, the relevant Feynman diagrams will be those of a 0^-theory. 

Moreover, the individual maximal interactions are identifiable by acting on special 
states: Assume that the Hamiltonian density maps a state 



There is a single interaction which achieves this particular reshuffling of spins: {2, 1, 3, 2}, 
cf. Fig. [TJ Therefore we could infer that c equals the coefflcient of this interaction, 
c = —12 — 2/^2,3 — 4e3a. The same is true for the other maximal interactions: If all lines 
going right are associated with cf) and the others with Z then a cj) will move past a Z 
towards the right at each elementary crossing. The effect will thus uniquely identify the 
corresponding interaction. 

Being the representation of the dilatation generator on the gauge invariant operators 
in the su(2) sector, the spin chain Hamiltonian in Tab. [T] is also the anomalous dimension 
matrix of operators in this sector. In any conformal field theory the eigenvalues of the 
anomalous dimension matrix are independent of the renormalization scheme. Its matrix 



ni\...(t)(i)Zz...) = c\...zz(j)(t)...) + ... . 



(2.7) 
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4a: {2,1,3,2} 



4b: {2,1,4,3} 



4c: {1,4,3,2} 



P 



P 




4d: {3,2,1,4} 



4e: {4,3,2,1} 



Figure 1: Maximal planar interactions up to four loops. Below the diagrams the permutation 
symbols are indicated. Solid and dashed lines correspond to two complex scalars in A/" = 4 
SYM. Above the diagrams suitable momenta to remove IR singularities are indicated. 

elements however do not generically have this property. At the level of the Hamiltonian 
in Tab. [H this is reflected by the fact that the undetermined coefficients e do not affect 
its eigenvalues [22]. 



3 Four-Loop Calculation 

With these preparations we are now in a position to compute the undetermined coef- 
ficients that appear in the spin chain Hamiltonian. Among the various approaches to 
finding the anomalous dimension matrix we shall consider the renormalization of compos- 
ite operators. It was successfully used in [27] to determine the two-loop and (under the 
assumption of proper BMN scaling for one-excitation BMN states) the all-loop dispersion 
relation. 

The renormalization of composite operators and the subtraction of sub divergences 
proceeds by introducing renormalization factors and counterterm diagrams analogous 
to the Bogoliubov R-operation. For our purpose this procedure was systematized in [3] 
where an iterative subtraction scheme was developed that allows the subtraction of entire 
subdiagrams. This is the scheme we shall use. 

We are therefore to compute Feynman diagrams with one vertex being the composite 
operator of interest and additional vertices dictated by the A/" = 4 SYM Lagrangian. As 
described in detail in the previous section, our goal is to find the entries of the four-loop 
anomalous dimension matrix that reshuffle scalar fields in a maximal way. Besides scalar 
fields, the internal lines of these diagrams may a priori also be fermions and gauge fields. 
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Two simple observations imply however that the situation is substantially simpler. 

By inspection of the diagrams in Fig. [1] it is easy to see that Feynman diagrams 
containing gauge fields cannot lead to such maximal reshuffling (in the sense described 
previously) of scalar fields. Indeed, using the fact that the gauge field interactions are 
fiavor-blind, one may see that replacing any of the four-point vertices by scalar-vector 
interactions leads to diagrams not exhibiting maximal reshuffling. 

Let us consider next scalar-fermion interactions. R-charge conservation implies that 
any diagram with external fermion fields and an insertion of an operator in the 5u(2) 
sector vanishes identically. Diagrams with internal fermion lines have a similar fate. 
To see this let us note that the Yukawa interactions of the A/" = 4 SYM Lagrangian 
are proportional to the SO (6) Dirac matrices. Since the fields of the su(2) sector are 
complex, the Dirac matrices appearing in these vertices will also carry complex vector 
indices. Their algebra, {F"', F^} = 7]°-^, implies that holomorphic matrices square to 
zero. Therefore the fiavor of scalar fields coupling to fermionic loops must alternate. It 
is then easy to see that for 2n fiavors at most n pairs can be interchanged. This does 
not lead to a maximal permutation and we can thus disregard fermion loops. 

The conclusion is therefore that the Feynman diagrams contributing to the entries 
of the four-loop anomalous dimension matrix describing a maximal reshuffling of spins 
are scalar diagrams in which each vertex contains two types of scalar fields and the 
interaction interchanges them. These are the diagrams listed in Fig. [TJ 

To compute the contribution to the anomalous dimension matrix we need to compute 
the amplitudes in Fig. [1] and isolate the overall ultraviolet divergence by subtracting all 
their UV sub divergences. While in general it is convenient to use a variant of dimensional 
regularization which preserves supersymmetry, in the context of our calculation making 
a definite choice is not an issue since all our diagrams have scalar internal lines. However, 
since all fields are massless we must be careful to separate the UV divergences from IR 
divergences. To this end we shall assign off-shell momenta to some of the external fieldsjl] 
they are chosen such that the number of momenta is minimal while all IR divergences 
are eliminated. It turns out that up to four-loop order it suffices that only two of the 
external fields carry momentum; an appropriate choice is depicted in Fig. [TJ 

All momentum integrals may be computed easily by reduction to a small set of master 
integrals. Common building blocks are bubble diagrams with arbitrary exponents for the 
two propagators; their expressions are (here and elsewhere the dimensionality of space- 



^The external fields can either belong to the operator being renormalized or be attached to the 
vertices of the A/" = 4 SYM Lagrangian. 

^We extract the overall momentum dependence for the sake of notational convenience. 



time is = 4 - 2e) I 




r(ai)r(a2) 



V{d - ai2) 
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where au = ai + 02. Indeed, diagrams 1, 2, 3b, 3c, 4d and 4e may be computed exactly 
by repeated identification of one-loop bubble subdiagrams. Once a bubble subintegral 
is evaluated, the exponent of the propagator carrying the momentum p flowing through 
the bubble is shifted by an integer multiple of the dimensional regulator. 

A similar iterative identification of bubble subintegrals reduces the diagrams 3a and 
4c to special cases of the two-loop master bubble integral 



T(ai,a2,a3,a4,a5) = {p ) 



2\ai2345— 



X 









as 







(3.2) 



with some arbitrary powers of propagators. For integer exponents ai_...^5 such integrals 
have been computed in the past (e.g. [2H])- We are however interested in situations 
when some of the exponents are not integers (special cases have been previously analyzed 
in [29]), being dependent on the dimensional regulator. Perhaps the most effective way of 
computing such integrals is to use the Mellin-Barnes parametrization [28j. The identity 



{a + by 



1 

27ri 



, T(-w)T(w 
FTT" 



(3.3) 



allows a straightforward evaluation of the Feynman parameter integrals and expresses 
the result of the momentum integral in terms of multiple contour integrals which can be 
evaluated through the residue theorem. This method has the advantage of producing 
explicit integral representations for the coefficients of the various powers of the dimen- 
sional regulator. The algorithm of [30j for the analytic continuation e — as well as the 
numerical evaluation of the resulting coeffcients has been successfully automated |31j . 
An MB parametrization of T(ai, 02, 03, 04, 05) is 



T(ai, 02, as, 04, as 



(47r) 



r(ai)r(a4)r(a5)r(rf - a^s 

dWidW2 r(ai45 - i + Wi2^ 



(3.4) 



(27ri)2 r(ai245 - f + Wu] 



T{-wi)T{-W2)T{a4 + Wi2)T{d - 01245 - ^^12) 



ai4 - wi)r(ai2345 -d + wi) r(| - a45 - W2)r(| -03 + ^2 



r( 



3d 
2 



ai2345 - Wi) 



T{a-i - W2) 



with the notation aijk... = ai + aj + ak + ... and similarly for Wijk...- While this parametriza- 
tion does not manifestly exhibit the symmetries of the original diagram, they are restored 
after the remaining integrals are performed. It is possible (though perhaps less effi- 
cient in terms of the necessary number of MB parameters) to construct a Mellin-Barnes 
parametrization manifestly exhibiting the (^2)^ symmetries of (13. 2p . The e-expansions 
of the two-loop integrals we shall require read: 

T(l, 1, 1, 1, e) = L(l, If (i + + K + (-i + TC(3))e^ + . . • 
T(l, 1, 1, 6, 1) = L(l, 1)^ (I + + fe' + (+f - f C(3))6^ + . 
T(l, 1, 1, 1, 26) = L(l, If (1 + |6 + le' + (-f + fm)e' + . . 
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T(l, 1,1,1 + e, e) = L(l, 1)^ + + + (+A + f C(3))e3 + . 

r(l, 1, 1, 2e, 1) = L(l, 1)2 (3^ + + §6^ + - f C(3))e=^ + ...). (3.5) 

Here it was convenient to factor out two powers of the one- loop bubble L(l, 1) which has 
the expansion 

L(l.l) = (4.)- ^'f?-^' (3,6) 

^ ' ^ ^ ' 167r2r(2-2e) ^ ' 



lek (1 + 2e + (4 - 3^.2)^2 + (8 - - |C(3))e^ + . . . 



In both of the last two integrals 4a and 4b it is trivial to isolate a factor L(l, 1). The 
remaining three-loop integrals may be computed in several ways. One approach makes 
use of integration by parts identities, known in this case as the triangle rule, to reduce 
them to combinations of one- and two-loop bubble integrals with various exponents (see 
App. ^for details). A second approach directly evaluates the three-loop integrals and 
in the process tests that the infrared region is non-singular. We list the necessary Mellin- 
Barnes integrals and the e-expansions of all diagrams in App. ^B] Needless to say, the 
two calculations lead to the same answer. 

For a vector of operators O, the relation between the bare and renormalized operators 
is given by the renormalization factor Z 

Qbare _ . Qren _ ^3 7^ 

The £-loop contribution to Z is found from the overall divergence of £-loop diagrams 
with exactly one insertion of a member of the vector O. To isolate the overall divergence 
it is necessary to include counterterm diagrams which are generated recursively by the 
lower-loop renormalization factor. These diagrams also eliminate the non-local momen- 
tum dependence. The relation between the renormalization factor and the anomalous 
dimension matrix (a.k.a. dilatation generator or spin chain Hamiltonian) Ti is standard: 

5n = \iuieZ-^—^ — Z. (3.8) 
e^o rfln^iYM 

This expression implies an exponential-like structure for the renormalization factor zjf] 
in particular, the derivative of Z must be left-proportional to Z and, in order that the STi 
be well defined in the e — *^ limit, the factor of proportionality can only have additional 
simple poles0 

A subtraction scheme that enforces these constraints and at each loop order isolates 
directly the contribution to the anomalous dimension matrix was described in for the 
use in two-point functions. An adapted version for the use in operator renormalization 
diagrams is presented in App. ^ where the explicit rules are given and then applied to 
the relevant diagrams. 



■^For operators which do not mix under RG flow the relation is Z = exp [e~^ dtt~^^{tgy^„)\ . 
^It is in principle possible that in a different renormalization scheme individual matrix elements could 
have divergent terms; however these terms should be removable by similarity transformations. 
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Finally, in this scheme the equation (13.81) reduces to the simple operation of picking 
the residue of the 1/e pole of the subtracted diagrams 



2(167r2)^ lime /(e) . (3.9) 



^0 



The factor of (IGvr^)^ corresponds to the normalization of the £-loop Hamiltonian and 
allows for a direct comparison of the quantity / to the coefficients in (12. ip . In our case 
this leads to 



h = 


-2, 


ha. = 


- 4 + 


4C(3), 


h = 


-2, 


hh = 


+10- 


12C(3), 


hsi = 


+4, 


he = 


+ 2 + 


8C(3), 


/3b = 


-4, 


hd = 


-10 + 


4C(3), 


4 = 


-4, 


he = 


-10, 





(3.10) 

which represent the coefficients relating the structures listed in Fig. [T] and the spin chain 
Hamiltonian (cf. Tab. [T]). Clearly, h, h, he and he reproduce the coefficients of {1}, 
{2, 1}, {3, 2, 1} and {4, 3, 2, 1}, respectively. The coefficients undetermined by symmetry 
considerations are fixed by our calculation to b^ 

^e2 = -l, esa = -2 - 3C(3), zes^ = -3 - C(3) (3.11) 

and 

P2,3 = 4C(3). (3.12) 

In particular we are able to uniquely fix the leading coefficient /32,3 for the dressing phase. 
It is in full agreement with the results of P ^ fT ^ ITTj. 



4 Conclusions and Outlook 

In this paper we have computed the four-loop dilatation operator in the su(2) sector of 
A/" = 4 SYM. The main observation which led to substantial technical simplifications is 
that the coefficients undetermined by symmetry constraints can be chosen to correspond 
to "maximal interactions" - i.e. interactions that reshuffie the spins in a maximal way. 
For appropriately chosen gauge theory operators these interactions are entirely deter- 
mined by Feynman diagrams with only scalar interactions. We found that, starting at 
four-loop order, the anomalous dimensions of long operators become transcendental; this 
may be traced to the dilatation operator acquiring transcendental coefficients. We have 
extracted the relevant coefficient of the dressing phase and found it identical to the one 
reproducing the four-loop cusp anomalous dimension computed in [l6l[T7]. Our result 
confirms the particular analytic continuation used to guess the dressing phase at weak 
coupling [T3] . 

^The conventional factors of i indicate that the Hamiltonian is not manifestly hermitian. With a 
proper choice of scalar product, however, it becomes quasi-hermitian as it should. 
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The main obstacle for computing higher-loop anomalous dimensions in any sector of 
A/" = 4 SYM and thus directly computing the S-matrix dressing phase is, as in all off-shell 
calculations, the proliferation of Feynman diagrams. In compact sectors the symmetries 
of the theory restrict (sometimes substantially) the structure of the anomalous dimension 
matrix. At any loop order the maximal interactions enjoy the same technical simplifi- 
cations as the ones employed in the calculations described here; moreover, it is possible 
that some of the relevant momentum integrals exhibit a recursive structure]^ It would be 
interesting to identify and compute them, thus providing a direct evaluation of important 
parts of the dressing phase. 
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A The Triangle Rule 

Consider the Feynman integral 



34 



F(ai, a2, as, a^, a^) = 




(A.l) 



This integral may be part of a larger Feynman diagram and the labels ai,...,5 represent the 
exponents of the propagators of the corresponding internal lines. Inserting the operator 
l^d/dl^ in the integral representing this diagram and equating the results of the action 
of the derivative on the original integrand and the result of the integration by parts leads 
to 

(d - as - as - 2a5)F = (as 2+ (5" - 1") + 3+ (5" - 4r)) F (A.2) 

where, for example, l^F(ai, as, as, a4, a^) = F{ai ± 1, as, as, a^, a^). 

The various terms in such a decomposition may acquire however spurious infrared 
divergences which are regularized by the dimensional regulator and - provided that the 
IR of the original integral was properly regularized - cancel when all terms are assembled. 

The triangle rule together with the straightforward evaluation of bubble integrals 



candidate for this property is { 



TO, . 



1, TO + 1, 



, 2, . . . , 71 + TO — 1, . . . , n — 1, n 



A- 
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leads to the following expressions for the diagrams in Fig. [T 



Il = 










h = 


ip'y 


^2^L(l,l)L(l + e,l), 




ha. = 


ip'y 


~3^L(l,l)r(l,l,l,l,e), 




/3b = 




-3^L(l,l)2L(l + e,l + e), 




/3c = 


ip'y 


-3^L(l,l)L(l + e,l)L(l + 2e, 1), 




ha = 


ip'y 


-4.^(1,1) 
l-3e 


L(2,e)T(l,l,l,2e,l)- 


L(2,3e)T(l,l,l,e,l)], 


hh = 


ip'y 


-^^L(l,l)2r(l,l,l,l + e,e), 




he = 


(p'y 


-4e^(l,l) 

1 -3e 


L(2, e)r(l, 1,1,1, 2e) - 


L(2,l)T(l,l,l,l + e,e 



hd 
he 



+ e L(l, 1 + e) T(l, 1,1,1, 2e) - e L(l, 1 + e) T(l, 1, 1, 2e, 1) 
(/)-^^ L(l, lfL{l + e, 1) L(l + 2e, 1 + e), 

(/)-^^ L(l, 1) L(l + e, 1) L(l + 2e, 1) L(l + 3e, 1). (A.3) 



Note that the integrals T(l,l,l,e, 1) and T(l, 1, 1, 2e, 1) can be evaluated further 
using the triangle rule 

T(l, 1, 1, e, 1) = (l(2, e) - L(2, 2e) + e L(l, 1 + e) - e L(l + e, 1 + e)) , (A.4) 



T(l,l,l,2e,l) = y^^(^L(2,2e) -L(2,3e) + 2eL(l,l + 2e) -2eL(l + e, l + 2e)). 



For compactness, we left them unexpanded in flA.3p : The integral ha can thus be evalu- 
ated as an analytic expression in e. 



B Three-Loop Integrals 

The three-loop integral that remains after one identifies a one-loop bubble subintegral 
in ha may be evaluated directly thus testing the application of the triangle rule and the 
correct infrared regularization of the contributions to the anomalous dimension matrix. 
A Mellin-Barnes parametrization of a master integral containing ha is 



/?M(ai,a2,a3) = (p')""^+'-=^"/' 




(47r 



-3d/2 



dwidw2dw^dw4 



X 



r(-wi)r(-w2)r(-w3)r(-w4)r(ai23 - j + W12) 

r(ai)r(a2)r(a3)r((i - ai23) 



^It is trivial to identify in the expressions of /4a and /4c the IR divergent components mentioned 
above. 
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r(f - ai3 - wi)r(| - ai2 - W2)T{ai + W12) 
^ r(l - u'i)r(3f - 3 - ai23 - W2)r{l + a^a - f + ^^^12) 
r(rf-2-ai23-^i23)r(f -l + u;3) 



X 



r(l-«;3) 



r((i - 2 - ai23 - W24)r(f + W24) 



r(3 + ai23 -d + W234)r(l + ai23 - 2 + «^1234)r(l - 2 ~ ^234) 
r(l - W3)T{d - 1 + ti;234) 

Then, 

h. = ipT^'Lil,l)BMil,e,l) (B.2) 

and its evaluation leads to the result hsted in equation (lB.6p . 
Similarly, he is a special case of the master integral 



5L(ai,a2,a3,a4,a5) = (^2)^12345 +3-3^/2 ^ 
A Mellin-Barnes parametrization is 




£"6 
3.2 a4 
as 




(B.3) 



5/7(01,02,03,04,05) = (47r) 



-3d/2 



dwidw2dw^dwi 



X 



r(-wi)r(-^2)r(i - f + 012 + ^i2)r(i + W12) 



r(oi)r(o2)r((i - 1 - 012) 

d 

2 



r(-«;3)r(-it;4)r(i - f + 034 + ^34)1(1 + w^^) 



r(o3)r(o4)r(rf- 1 -034) 
X r(f - 1 - oi - wi)r(f - 1 - 02 - W2)r(f - 1 - 03 - W3)r(f - 1 - 04 - ^4) 
^ r(2 + 0123456 - f + W13) r(f - 05 + W24) r(f - 2 - 012345 - ^1234) 

T{2d - 2 - 0123456 - Wi^) r(05 - W24) r(2 + 012346 -d + 1^1234) 

Then 

he = (p')"^' L{1, 1) lim BL{1, e, 1, 1, i/, 1) (B.5) 

whose evaluation leads to the result listed in equation (1B.6p . 

The integrals listed here are useful for the calculation of the coefficients of the higher- 
loop Hamiltonian in the su(2) sector. The resulting e-expansions of the integrals in Fig. [T] 
read 

^1 = 1^ (^) ^ (l + 2e + (4 - ^vr^)e^ + (8 - ^ - 1^)^ + ...), 
h = {h?{\ + \e + \e'+{l-?>a^)y + --] 

^3a = (/if(| + |^ + F+(-| + ¥C(3)y + 
/3b = (/if (I + f 6 + P + (f - f C(3))e=' + . . 
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^3c=(/if(| + ie + fe^+(f-^C(3))e^ + ...; 

/4b=(/ir(^ + ^e+ie^+(A + fC(3))e^ 
/4c=(/ir(| + ie + p+(f-|C(3))e3 + ...; 

/4d=(/ir(| + ie + fe^+(f-fC(3))e^ + 
he = + + fe^ + - f C(3))6^ + ...). (B.6) 

C Subtraction Scheme 

Here we describe the subtraction scheme used to extract the contributions to the anoma- 
lous dimensions without having to insert couterterms at each stage of the calculation. 
For each connected diagram drawn with the composite operator as the lower-most vertex, 
one 

• partitions it in all possible connected subdiagrams (including the trivial partition into 
a single subdiagram) and interprets those diagrams as contributing to the renormal- 
ization of a composite operator, 

• discards all partitions which are interconnected horizontally (all partial diagrams 
must be "dropped" onto the composite operator from above in a well-defined se- 
quence, in similarity to a famous arcade game), 

• discards all partitions for which there are two or more top-most diagrams, 

• evaluates the momentum integrals for the of the remaining partitions, 

• sums the products of the momentum integrals for each partition weighted by £(—!)"■ 
{n is the number of partial diagrams and i is the loop number of the top-most diagram 
in the partition) J§ 

Applying this scheme to the diagrams in Fig. [T] we find the following subtracted 
integrals / 

h = -/i, 

h = -2/2 + II 

h. = -3/3a + 2/2/1, 

/sb = -3/3b + 4J2/1 - II 
he = "~3/3c + 3/2/1 /i ) 

/4a = -4/4a + hJl + ^hhh + " ^hll + /' , 

/4b = -4/4b + 3/3a/l + /Sb/l + 2ll - 2hll 



®It is easy to see that this weight is related to a derivative with respect to the loop-counting parameter, 
as in equation 
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he = -4/4c + hJl + hch + 21^ - hll 

Im = -4/4d + 3/3b/i + 3/3c/i + 2/2 - 5/2/1' + It, 

he = -4/4e + 4/3c/l + 2/| - 4/2/1' + A'- (CI) 

It is not hard to find that the the quantities / exhibit only a simple pole in the e- 
expansion. A strong crosscheck of the correctness of the subtraction is the cancellation 
of non-local and divergent momentum- dependence that arises in (IB. 61) . 
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